Rotating Boson Stars and Q-Balls II: 
Negative Parity and Ergoregions 
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We construct axially symmetric, rotating boson stars with positive and negative parity. Their flat 
space limits represent spinning Q-balls. Q-balls and boson stars exist only in a limited frequency 
range. The coupling to gravity gives rise to a spiral-like frequency dependence of the mass and 
charge of boson stars. We analyze the properties of these solutions. In particular, we discuss the 
presence of ergoregions in boson stars, and determine their domains of existence. 
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Q-balls represent stationary localized solutions in flat space. They arise, when a complex 
scalar field has a suitable self- interaction [1, 2]. The global phase invariance of the scalar 
field theory is associated with a conserved charge Q, corresponding for instance to particle 
number [1]. 

The time-dependence of the Q-ball solutions resides in the phase of the scalar field and 
is associated with a frequency tu s . Q-balls exist only in a certain frequency range, w m i n < 
ui s < w max , determined by the properties of the potential [1-5]. At a critical value of the 
frequency, both mass and charge of the Q-balls assume their minimal value, from where they 
rise monotonically towards both limiting values of the frequency. Considering the mass of 
the Q-balls as a function of the charge, there are thus two branches of Q-balls, merging and 
ending at the minimal charge and mass. 

The simplest type of Q-balls is spherically symmetric. These possess finite mass and 
charge, but carry no angular momentum. Besides the fundamental spherically symmetric 
Q-balls there are also radially excited Q-balls [4]. The fundamental spherically symmetric 
Q-balls are stable along their lower mass branch, as long as their mass is smaller than the 
mass of Q free bosons [1]. 

Recently, the existence of rotating Q-balls was demonstrated [4, 5]. The rotation is 
achieved by including an additional dependence of the phase factor of the scalar field on 
the azimuthal angle (p, where the proportionality constant n must be integer. The resulting 
stationary localized Q-ball solutions then possess finite mass and finite angular momentum. 
Interestingly, their angular momentum J is quantized in terms of their charge Q, J = nQ 
[4, 6]. Rotating Q-balls with n = 1 thus have the smallest angular momentum J for a given 
charge Q. The energy density and angular momentum density of rotating Q-balls possess 
axial symmetry. There are no (infinitesimally) slowly rotating Q-balls [4] . 

The field equations allow for Q-balls with positive and negative parity. The scalar field 
is then symmetric resp. antisymmetric w.r.t. reflections [4]. For a given charge Q one finds 
thus two sequences of rotating Q-balls, whose angular momentum increases with n: the 
positive parity sequence n + and the negative parity sequence n~ . While the energy density 
of rotating Q-balls with positive parity corresponds to a torus, the energy density of rotating 
Q-balls with negative parity corresponds to a double torus. 

When the scalar field is coupled to gravity, boson stars arise [3, 7-9]. The presence 
of gravity has a crucial influence on the domain of existence of the classical solutions. 
Stationary spherically symmetric boson stars (0 + ) also exist only in a limited frequency 
range, ujq(k) < uj s < w max , where n denotes the strength of the gravitational coupling. They 
show, however, a rather different type of frequency dependence. In particular, in the lower 
frequency range, the boson star solutions are not uniquely determined by the frequency. 
Instead a spiral-like frequency dependence of the charge and the mass is observed, where 
the charge and mass approach finite limiting values at the center of the spiral. Furthermore, 
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the charge and the mass of these boson stars tend to zero, when the maximal value of the 
frequency is approached [3, 7]. 

Rotating boson stars with positive parity have been obtained before [5, 6, 10-12]. A study 
of the frequency dependence of the 1 + boson stars showed, that their frequency dependence 
is analogous to the one of the non- rotating + boson stars [5]. In particular, their charge 
and mass also exhibit a spiral structure. Like rotating Q-balls, rotating boson stars exhibit 
the angular momentum quantization, J = nQ [6]. 

Here we reconsider rotating boson stars with positive parity and, for the first time, present 
rotating boson stars with negative parity. We construct families of boson stars of both 
parities numerically for constant values of the gravitational coupling strength and n = 1 
and 2. We investigate the properties of these solutions. In particular, we analyze their 
frequency dependence, and we demonstrate the appearence and development of crgorcgions 
for these families of boson star solutions [13]. 

In section II we recall the action, the general equations of motion and the global charges. 
In section III we present the stationary axially symmetric Ansatz for the metric and the 
scalar field, we evaluate the global charges within this Ansatz, and present the boundary 
conditions for the metric and scalar field function. We discuss Q-ball solutions and their 
properties in section IV and boson star solutions and their properties in section V. Section 
VI gives our conclusions. 

II. ACTION, EQUATIONS AND GLOBAL CHARGES 
A. Action 



We consider the action of a self-interacting complex scalar field $ coupled to Einstein 
gravity 



V=gd 4 x , (l) 



where R is the curvature scalar, G is Newton's constant, the asterisk denotes complex 
conjugation, 

and U denotes the potential 

U{\*\) = A|$| 2 (|$| 4 - a|$| 2 + b) = A (0 6 - a<f + btf) , (3) 

with |$| = (f>. The potential is chosen such that nontopological soliton solutions [1], also 
referred to as Q-balls [2], exist in the absence of gravity. As seen in Fig. 1, the self-interaction 
of the scalar field has an attractive component, and the potential has a minimum, U(0) = 0, 
at $ = and a second minimum at some finite value of |$|. The boson mass is thus given by 
tub = \/A6. In the numerical calculations of the solutions presented, we choose the potential 
parameters [4, 5] 

A=l, a = 2, 6=1.1. (4) 



B. Equations 

Variation of the action with respect to the metric leads to the Einstein equations 

G>„ = R^v — —g^vR = kT^ u , (5) 
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FIG. 1: Left: The potential U(<j>) versus <f> for A = 1, a = 2 and b = 1.1 resp. 6=1. Right: The 
effective potential V(4>) = \ oj 2 <f> 2 — \U(<t>) versus <f> for several values of the frequency w s . 



with k = 8irG and stress-energy tensor T^j, 
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Variation with respect to the scalar field leads to the matter equation 

dU 



□ + 



$ = o 



(6) 
(7) 

(8) 



where □ represents the covariant d'Alembert operator. Equations (5) and (8) represent the 
general set of non-linear Einstein-Klein-Gordon equations. 



C. Global Charges 

The mass M and the angular momentum J of stationary asymptotically fiat space-times 
can be obtained from their respective Komar expressions [14], 



and 



M = — / R^n^dV , (9) 



E 



J = ~8^G I n '"'""'>"' ,V ' 



e 



Here S denotes an asymptotically flat spacelike hypersurface, n M is normal to S with n^n^ = 
— 1, dV is the natural volume element on S, £ denotes an asymptotically timclikc Killing 
vector field and 77 an asymptotically spacelike Killing vector field [14]. Replacing the Ricci 
tensor via the Einstein equations by the stress-energy tensor yields 

M = 2 jf (t^ - l - 9lMV Tj) n»?dV , (11) 



and 



J = - I ( '/;„■ ^ '/„■ /•' )n^dV . (12) 
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A conserved charge Q is associated with the complex scalar field $, since the Lagrange 
density is invariant under the global phase transformation 

$ -> $e M , (13) 

leading to the conserved current 

f = -* - $0"$*) , j% = 0. (14) 

III. ANSATZ AND BOUNDARY CONDITIONS 
A. Ansatz 

To obtain stationary axially symmetric solutions, we impose on the space-time the pres- 
ence of two commuting Killing vector fields, £ and rj, where 

£ = v = d v (15) 

in a system of adapted coordinates {t, r, 9, ip}. In these coordinates the metric is independent 
of t and <p, and can be expressed in isotropic coordinates in the Lewis-Papapetrou form [15] 



ds 2 = -f dt2 + J 



h (dr 2 + r 2 d9 2 ) + r 2 sin 2 6 (dtp - ^ dtj ' 



(16) 



The four metric functions f, I, h and u> are functions of the variables r and only. 

The symmetry axis of the spacetime, where r\ — 0, corresponds to the z-axis. The 
elementary flatness condition [16] 

y Y' A* 

^^ = 1, X = ^» (17) 

then imposes on the symmetry axis the condition [17] 

h\ e=0 = h\ e=n = 1 . (18) 

For the scalar field $ we adopt the stationary Ansatz [6] 

$(t,r,6,<p) = <j>(r,6)e i0J ° t+in f , (19) 

where (f>(r, 9) is a real function, and uj 8 and n are real constants. Single-valuedness of the 
scalar field requires 

= $(2tt + ip) , (20) 

thus the constant n must be an integer, i.e., n = 0, ±1, ±2, .... We refer to n as the 
rotational quantum number. When n ^ 0, the phase factor exp (imp) prevents spherical 
symmetry of the scalar field <f>. 

Solutions with positive and negative parity satisfy, respectively, 

0(r,7r-0) - cj>(r,6) (21) 
0(r,7r-0) - -0(r,0) , (22) 

To construct stationary axially symmetric boson star solutions a system of five coupled 
partial differential equations needs to be solved. In contrast, for Q-balls the metric is the 
Minkowski metric, i.e.,/ = l = h = l,u = 0. Here, at least in principle, only a single partial 
differential equation for the scalar field function needs to be solved. 
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B. Mass, angular momentum and charge 



The mass M and the angular momentum J can be read off the asymptotic expansion of 
the metric functions / and ui, respectively, [15] 



2JG „ 



i.e., 



M = lim r 2 d r / , J = lim r 2 u . 

ZG r^oo ZG r^oc 



(23) 



(24) 



This is seen by considering the Komar expressions Eqs. (9) and (10), with unit vector 
= (1,0,0, u/r)/VJ, and volume element dV = l/VT \g\ ^ 2 dr d6 dip, leading to [18] 



M 



= -g^ J I%V=gdrd0dip 



lim — 

r^oo 8irG 
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f 



r 2 sin ( 



and similarly 



J = lim 



df I . 2fl 

sin' 9lu 

dr f 

I 3 / 2 



duj lo 
dr r 



dO 



r 2 sin 3 I 



dto 
dr 



d9 . 



(25) 



(26) 



Insertion of the asymptotic expansions of the metric functions then yields expressions (24) . 

Alternatively, the mass M and the angular momentum J can be obtained by direct 
integration of the expressions (11) and (12), where 

= J (27* -T£) Ig^drdOdcp, (27) 
corresponds to the Tolman mass, and 

J = - J T^g^drdOckp. (28) 
The conserved scalar charge Q is obtained from the time-component of the current, 
Q = - j f \g\ 1/2 drdOdip 



= Atxuj s 



OO P7T 



|l/2_ 



1 + 



drdO 



From Eq. (28) for the angular momentum J and Eq. (29) for the scalar charge Q, one 
obtains the important quantization relation for the angular momentum, 



J = nQ , 



(30) 



first derived by Schunck and Mielke [6], by taking into account that = nj*, since <9 V $ 
m$. Thus a boson star with n = carries no angular momentum, J = 0. 



C. Boundary conditions 



The choice of appropriate boundary conditions must guarantee that the boson star solu- 
tions are globally regular and asymptotically flat, and that they possesses finite energy and 
finite energy density. 
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For rotating axially symmetric boson stars appropriate boundary conditions must be 
specified for the metric functions f(r,9), l(r,6), h(r,9) and u(r,6), and the scalar field 
function <f>(r,9) at infinity, at the origin, on the positive z-axis (9 = 0), and, exploiting the 
reflection properties w.r.t. 9 — > n — 9, in the xy-plane (9 = n/2). 

For r — > oo the metric approaches the Minkowski metric 77^,3 and the scalar field assumes 
its vacuum value <E> = 0. Accordingly, we impose at infinity the boundary conditions 

At the origin we require 

d r f\ r=0 = , d r l\ r=0 = , h\ r =o = 1 , uj\ r =o = , (f)\ r =0 = . (32) 

Note, that for spherically symmetric boson stars the scalar field has a finite value 0o at the 
origin, though. 

For 9 = and 9 = ir/2, respectively, we require the boundary conditions 

d e f\e=o = , d e l\ g =o = , h\ g=0 = 1 , d e uj\e=o = , </)\e=o = , (33) 

and 

def\e=7r/2 = , d e l\e^/ 2 = , d e h\ g=7T / 2 = , d e u\ g=7T /2 = , d e <p\ g=7T /2 = , (34) 

for even parity solutions, while for odd parity solutions (f>\e=7r/2 = 0. 

In non-rotating solutions, the metric function ut is identically zero. In spherically symmet- 
ric solutions the functions depend on the radial coordinate only, with metric function h = 1 . 
For the flat space solutions the metric corresponds to the Minkowski metric everywhere. 



D. Numerical method 



Rotating solutions are obtained when n^O. The resulting set of coupled non-linear par- 
tial differential equations is solved numerically [19] subject to the above boundary conditions, 
Eqs. (31)-(34). Because of the power law fall-off of the metric functions, we compactify space 
by introducing the compactified radial coordinate 

r = -!— . (35) 

The numerical calculations are based on the Newton-Raphson method. The equations are 
discretized on a non-equidistant grid in f and 9. Typical grids used have sizes 90 x 70, 
covering the integration region < f < 1 and < 9 < it/2. 



IV. Q-BALLS 

A. Spherically symmetric Q-balls 

We first briefly review the main features of spherically symmetric Q-balls (see e.g. [1, 2, 4], 
or [3] for a review). The equation of motion for spherically symmetric Q-balls is given by 

It may be interpreted to effectively describe a particle moving with friction in the effective 
potential 

V(cP) = ^W- 1 -U(4>), (37) 

exhibited in Fig. 1. 
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Q-balls exist in the frequency range 

^min <"s< W L . (38) 

where the maximal frequency oj max and the minimal frequency w m i n are determined by the 
properties of the potential U [4] , 

u 2 s < ^ ax = l -U"{0) = \b, (39) 

lu 2 > u£ in = mm [t/(0)/</> 2 ] = \ (b-^j . (40) 
The mass M of the Q-balls can be obtained from 

/'OC 

M = 4it [lu 2 s 4> 2 + <t> 12 + U((f>)]r 2 dr , (41) 

where the prime denotes differentiation with respect to r, and and their charge Q from 

Q(u>s) = 8iru s cj) 2 r 2 dr . (42) 
Jo 

Fixing the value of uu s in the allowed range, one then obtains a sequence of Q-ball solutions, 
consisting of the fundamental Q-ball and its radial excitations [4] . The boson function <f> of 
the fundamental Q-ball has no nodes, while it has k nodes for the fc-th radial excitation. 

We exhibit the fundamental Q-ball solution (k = 0) and the first radial excitation (k = 1) 
in Fig. 2. Here the charge Q is shown versus the frequency lu s . As seen in the figure, at a 
critical value w cr , which depends on k, the charge assumes its respective minimal value Q cr . 
Thus two branches of solutions exist for each k, the lower branch u s < iv a and the upper 
branch uj s > w cr , which merge and end at w cr . When uj s — ► w min as well as when uo s — » o> max 
the charge diverges [1]. 

In Fig. 2 also the mass M of these two sets of solutions is shown, but it is now exhib- 
ited versus the charge Q. The lower mass branches correspond to the lower values of the 
frequency, uo s < oj ct , while the upper mass branches correspond to higher values of the fre- 
quency, u> s > o-> C r- As long as the mass is smaller than the mass of Q free bosons, M < tubQ, 
the fundamental Q-ball solutions are stable, both classically and quantum mechanically [1]. 
Classically the solutions remain stable all along the lower branch. The classical stability 
changes only at the critical point, Q CT , where the solutions acquire an unstable mode [1]. In 
contrast, from a quantum perspective the stability changes when the lower mass branch in- 
tersects the line corresponding to the mass of Q free bosons. On the unstable upper branch 
the mass approaches the mass M = iubQ of Q free bosons from above, as uj s — + w max [1]. 

We remark, that also angularly excited non-rotating Q-balls have been obtained recently 
[20]. These new solutions are axially symmetric, and their angular dependence is given by 
spherical harmonics. They have higher energy and charge than the spherically symmetric 
Q-balls, but lower energy and charge than the rotating Q-balls. In addition, interacting 
Q-balls have been studied [20]. 



B. Rotating Q-balls 

The existence of rotating Q-balls was shown by Volkov and Wohnert [4] . Based on the 
Ansatz (19) for the scalar field $ [6], rotating Q-balls are solutions of the field equation 



d 2 , 2 d t 1 d 2 cos9 d n 2 2 \ x _ ldU(<j>) 



dr 2 rdr r 2 d6 2 r 2 sm6d6 r 2 sin 2 6 s ) 2 
with mass [4] 



(43) 



M = 2ir ^ P (lu 2 ^ 2 + {d^f + \(dg^) 2 + f^ 2 - + U(ct>)) r 2 dr sm6 d6 , (44) 
Jo Jo V r r 2 sm z 6 J 
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FIG. 2: Left: The charge Q versus the frequency u s for spherically symmetric fundamental Q-balls 
(k = 0) and their first radial excitations (k — 1). Also shown are the limiting values of the frequency, 
Wmin and Wmax. Right: The mass M versus the charge Q for spherically symmetric fundamental 
Q-balls (k = 0) and their first radial excitations (k — 1). Also shown is the mass for Q free bosons, 
M = tubQ- The mass on the upper branches is very close this mass. The k — branches close to 
the critical charge Q cr are seen in the inset. 



and charge 

/>OC />7T 

Q = 4ttuj s / 4> 2 r 2 dr sm.6d6 . (45) 
Jo Jo 

Their angular momentum satisfies the quantization relation J = nQ. 

In their pioneering study [4] Volkov and Wohnert showed, that for a given value of n 
there are two types of solutions, positive parity Q-balls n + and negative parity Q-balls n~ . 
In the following we restrict our considerations to fundamental rotating Q-balls. 



1. Positive parity Q-balls 

We first recall the properties of rotating Q-balls with positive parity and, at the same 
time, extend the previous studies [4, 5] to higher rotational quantum number n. In Fig. 3 
we illustrate the scalar field 4> and the energy density T u of rotating Q-balls with charge 
Q = 410, n p = 1 + , 2 + and 3 + . For rotating Q-balls the scalar field <p vanishes on the 
2:-axis. The energy density T u of even parity rotating Q-balls typically exhibits a torus- 
like structure. This torus-like shape becomes apparent by considering surfaces of constant 
energy density of e.g. half the respective maximal value of the energy density. As seen in 
the figure, the maximum of the energy density increases and shifts towards larger values of 
the coordinate p = rsm9, when n and thus the angular momentum increases. In terms of 
the energy density tori this manifests as an increase of their radii in the equatorial plane 
with increasing n. (Note, that in [5] a weighted energy density was shown.) 

We next address the frequency dependence for these rotating Q-ball solutions with positive 
parity. In Fig. 4 we exhibit the charge Q versus the frequency lo s for Q-balls with n = — 3. 
We observe the same upper limiting value w max , Eq. (39), for the frequency u s , as for non- 
rotating Q-balls, ensuring asymptotically an exponential fall-off of the scalar field (p. For 
a given frequency uj s the charge of a Q-ball increases with n. We thus conclude, that the 
frequency of rotating Q-balls is limited by the minimal frequency w m i n , Eq. (40), as well, 
independent of n. 

The mass of these rotating Q-balls with positive parity shows the same critical behaviour 
and thus cusp structure as the mass of the non-rotating Q-balls. This is also illustrated 
in Fig. 4. Thus we conclude, that these rotating solutions are classically stable along their 
lower mass branches. Overall we note, that the sets of rotating Q-balls with positive parity 
exhibit the same general pattern as the sets of non-rotating Q-balls. 
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FIG. 3: The scalar field <f> (left column) and energy density T t t (right column) for rotating positive 
parity Q-balls with charge Q = 410 and n p — 1 + (upper row), 2 + (middle row) and 3 + (lower row) 
versus the coordinates p = r sin 6 and z — r cos 9 for z > 0. 



2. Negative parity Q-balls 

Turning to rotating Q-balls with negative parity, we first illustrate in Fig. 5 the scalar 
field 4> and the energy density T u of Q-balls with fixed charge Q = 720 and n p = 1~ and 
2~ . For these Q-balls the scalar field <fr vanishes not only on the z-axis but everywhere 
in the xy-plane, as well. The energy density T u now typically exhibits a double torus-like 
structure, where the two tori are located symmetrically w.r.t. the equatorial plane [4]. When 
n and thus the angular momentum increases, the maxima of the energy density again shift 
towards larger values of the coordinate p — rsmO. For the double tori of the energy density 
this shift with increasing n manifests as an increase of the radii of the tori in the equatorial 
plane. 

In Fig. 6 we exhibit the charge Q versus the frequency lo s for negative parity Q-balls with 
n p = 1~ and 2~, and compare to the positive parity Q-balls. While we conclude, that 
for the negative parity Q-balls the limiting values ui m i n and ui max are retained, numerical 
difficulties have unfortunately prevented a more complete study. 

Comparing the sets of negative and positive parity Q-balls, we conclude, that they exhibit 
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FIG. 4: Left: The charge Q versus the frequency w s for Q-balls with n p = + , 1 + , 2 + , 3 + . The 
asterisks mark the critical values of the charge. Right: The mass M versus the charge Q for Q-balls 
with n p = 0+, 1+, 2+, 3+. The upper branches of the mass M are hardly discernible (on this scale) 
from the mass of Q free bosons, M = msQ- 





FIG. 5: The scalar field <fr (left column) and energy density T tt (right column) for rotating negative 
parity Q-balls with charge Q = 720 and n p = 1~ (upper row) and 2~ (lower row) versus the 
coordinates p — r sin 6 and z = r cos 9 for z > 0. 



the same general pattern. However, concerning the magnitudes of their charges and masses, 
the positive parity Q _ balls and the negative parity Q-balls do not quite alternate. Instead, 
they are ordered according to n p = + , 1 + , 2+, 1~, 3 + , 2~, as seen in Fig. 6. Thus for 
Q-balls of a given charge Q, the double torus- like structure of the negative parity solutions 
leads to a considerably higher mass than the single torus-like structure of the positive parity 
solutions. 
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FIG. 6: Left: The charge Q versus the frequency w s for negative parity Q-balls with n p = 1~ and 
2~. For comparison the positive parity Q-balls with n p — + , 1 + , 2 + , 3 + are also shown. Right: 
The mass M versus the charge Q for the same set of solutions together with the mass of Q free 
bosons, M — msQ- 



V. BOSON STARS 



A. Spherically symmetric boson stars 



When the scalar field is coupled to gravity, boson stars arise (see e.g. [8, 10] for reviews). 
We first demonstrate the effects of gravity for spherically symmetric boson stars in Fig. 7, 
where we exhibit the charge Q versus the frequency ui s for several values of the gravitational 
coupling constant n. 
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FIG. 7: Left: The charge Q versus the frequency lu s for fundamental spherically symmetric boson 
stars (k = 0) in the full range of existence for a sequence of values of the gravitational coupling 
constant k. Right: The charge Q versus the frequency uj s in the frequency range of the spiral for 
the gravitational coupling constant k = 0.2. 



For solutions in curved space the frequency uj s is also bounded from above by w max , 
Eq. (39), ensuring an asymptotically exponential fall-off of the scalar field. But, unlike for 
flat space solutions, the charge Q does not diverge for boson stars, when oj s approaches 
w m ax- Instead the charge tends to zero in the limit uj s — > uj max . 

Also, for the smaller values of ui s the solutions approach the limiting lower value w m i n , 
Eq. (40), no longer monotonically. Instead a spiral-like behaviour arises in the presence of 
gravity, where the boson star solutions exhibit an inspiralling towards a limiting solution 
[7]. The location and size of the spiral depend on the gravitational coupling strength k, and 
can be well below w m j n for small values of n (see Fig. 7) [5] . 

The limit n — > oo is obtained by introducing the scaled scalar field <j>(r) = y/n(p(r) [5]. 
Substituting <fi(r) in the field equation and taking the limit k — > oo, the terms non-linear 
in 4>{r) vanish. Similarly, in the Einstein equations the terms of higher than second order 
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FIG. 8: Left: The scaled charge nQ versus the frequency u> 3 for fundamental spherically symmetric 
boson stars in the limit k — ► oo. For comparison the scaled charge kQ is also shown for k = 0.2. 
Right: The mass M versus the charge Q for fundamental spherically symmetric boson stars for the 
value of the gravitational coupling constant k = 0.2 

in (f>(r) vanish and the dependence on k cancels. The resulting set of differential equations 
is then identical to the original one, except that k = 1 and U((f>) = \b(f> 2 7 reducing the 
potential to a mass term. This also yields a spiral pattern for the (scaled) charge Q = kQ, 
as seen in Fig. 8. 

The mass M has an analogous dependence on the frequency lu s as the charge Q. When 
the mass M is considered as a function of the charge Q, however, one now observes a new 
intricate cusp structure [5, 7]. As gravity is weekly coupled, the single cusp present in 
flat space (associated with the minimal value Q m [ n of the charge) is retained. But it is 
supplemented by a second cusp, since in the presence of gravity the charge and the mass 
do not diverge for uj s — ► w max , but tend to zero in this limit. For larger values of the 
gravitational coupling k these two cusps merge and disappear. However, an additional set 
of cusps arises in the presence of gravity, which is due to the appearance of the spirals, and 
therefore only present in curved space [5, 7]. The cusp structure is illustrated in Fig. 8. 

The classical stability of boson stars can be analyzed according to catastrophe theory, 
implying a change of classical stability at each cusp encountered [21]. In the following we 
refer to the main branch of solutions as the branch between the maximum of the charge 
and the mass (at a relatively small value of the frequency u> 3 ) and the local minimum of 
the charge and the mass (at a relatively large value of the frequency u s ). Note, that for 
the flat space solutions the minimum is a global minimum, the charge and the mass assume 
their critical values here, while for large gravitational coupling k the local minimum has 
disappeared, and the relevant minimum corresponds to the global minimum at w max . 

Along this main branch the boson star solutions are classically stable. We note though, 
that a small part of this main branch may become quantum mechanically unstable, when 
M > iubQ close to the first cusp. Since at each cusp terminating the main branch a negative 
mode is acquired, the solutions on the adjacent branches become classically unstable. In the 
spiral, the boson stars are expected to acquire at each further cusp an additional negative 
mode, making the solutions increasingly unstable. The physically most relevant set of 
solutions should thus correspond to the stable boson star solutions on the main branch. 

B. Rotating Boson Stars 

We now turn to rotating boson stars, which emerge from rotating Q-balls, when the 
gravitational coupling constant is increased from zero. We first recall the previous results 
on rotating positive parity boson stars with n p = 1 + [5, 6, 12] and present new results 
for n p = 2+. We then turn to rotating negative parity boson stars. For both cases we 
demonstrate and analyze the occurrence of ergoregions. We here restrict our analysis to 
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fundamental rotating boson stars. 



1. Positive parity boson stars 



We begin the discussion by recalling the main results on boson stars with n p = 1 + , 
obtained previously [5]. For that purpose we exhibit in Fig. 9 the charge Q of 1 + boson 
stars versus the frequency uj s for several values of the gravitational coupling constant re, 
including the flat space limit. 
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FIG. 9: Left: The charge Q versus the frequency u) s for rotating 1 + boson stars for several values 
of the gravitational coupling constant re. Also shown are the limiting flat space values. The dots 
indicate the onset of ergoregions. Right: The mass M versus the frequency lu 8 for the same set of 
solutions. 



Again, the frequency u) s of the solutions is bounded from above by w max , Eq. (39), ensuring 
an asymptotically exponential fall-off of the scalar field. Also, the charge of the rotating 
boson stars tends to zero, when w max is approached. Likewise, 1 + boson stars exist only in 
a finite interval of the frequency oj s , where the lower bound increases with increasing re. 

Furthermore, we observe the presence of spirals also for these rotating boson stars, i.e., the 
1 + boson stars exhibit an inspiralling towards a limiting solution just like the non-rotating 
boson stars. For small values of the coupling constant re the spirals are located in the region 
of small frequency uj s , for larger coupling re the spirals shift towards larger frequencies. 
A good and rather complete determination of the spirals is, however, numerically often 
extremely difficult [5]. 

The mass M of the 1 + boson stars has an analogous dependence on the frequency oj s 
as the charge Q, as seen in Fig. 11. When the mass M of these rotating boson stars is 
considered as a function of the charge Q, we observe an analogous cusp structure as for 
the non-rotating boson stars [5]. In particular, for a given gravitational coupling re, a main 
branch of rotating boson star solutions is present, and located between the maximum of 
the charge and the mass (at a relatively small value of the frequency u> s ) and the (for small 
gravitational coupling only local) minimum of the charge and the mass (which is not part 
of the spiral). Following arguments from catastrophe theory, the boson star solutions are 
again expected to be classically stable along this main branch. 

Let us now turn to positive parity boson stars with higher quantum number n and thus 
greater angular momentum. We exhibit in Fig. 10 the charge Q and the mass M of 2 + bo- 
son stars versus the frequency ui s for several values of the gravitational coupling constant re, 
including the flat space limit. Because these calculations are numerically very involved and 
time-consuming, we have concentrated on the main branches of these 2+ boson stars, omit- 
ting the numerically even more involved spirals. Clearly, the 2 + solutions show a completely 
analogous re- and frequency dependence (in the regions studied) as the 1 + solutions. 

To further demonstrate the dependence of the positive parity boson stars on the rotational 
quantum number n, we exhibit in Fig. 11 the charge Q and the mass M as functions of 
the frequency uj s for boson stars with quantum numbers n p = + , 1 + , 2 + and gravitational 
coupling re = 0.2. 
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FIG. 10: Left: The charge Q versus the frequency ui s for rotating 2 + boson stars for several values 
of the gravitational coupling constant k. Also shown are the limiting flat space values. The dots 
indicate the onset of ergoregions. Right: The mass M versus the frequency lo s for the same set of 
solutions. 
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FIG. 11: The charge Q (left column) and the mass M (right column) versus the frequency ui s for 
boson stars with n p — + , 1 + , 2 + for the gravitational coupling k = 0.2. 



Considering the spatial structure of the solutions, the scalar field cf) typically gives rise 
to a torus-like energy density T tt for the rotating 1+ and 2+ boson stars, just as for the 
rotating n + Q-balls. Again, the torus-like shape becomes apparent by considering surfaces 
of constant energy density of e.g. half the respective maximal value of the energy density. As 
in the flat space solutions, the maximum of the energy density of the boson stars increases 
and shifts towards larger values of the coordinate p = rsin0, when n and thus the angular 
momentum increases. Thus in terms of the energy density tori this again manifests as an 
increase of their radii in the equatorial plane with increasing n. 

Rotating objects may possess ergoregions. But unlike the ergoregions of black holes, the 
ergoregions of regular objects like boson stars would signal the presence of an instability 
[13, 22-24]. This instability can be traced back to superradiant scattering, and its possible 
presence was put forward by Cardoso et al. [13] as an argument against various black hole 
doubles as potential horizonless candidates for compact dark astrophysical objects. 

Therefore we now analyze the occurrence of ergoregions for rotating positive parity boson 
stars. The ergosurface is defined by the condition 



fit i = 



-f+j sin 2 



(46) 



in the metric parametrization Eq. (16). The ergoregion resides inside this ergosurface. The 
presence of an ergosurface thus implies instability of the respective boson stars [13]. 

Examining the condition g tt = for the sets of boson star solutions, we indeed observe 
that rotating boson stars possess ergoregions in a large part of their domain of existence. 
To discuss these ergoregions, let us recall Figs. 9 and 10, where the charge Q is exhibited 
versus the frequency ui s for several sets of boson stars solutions with n p = 1 + and 2 + , 
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respectively. The dots in these figures indicate the onset of ergoregions and thus ergoregion 
related instability at a corresponding critical frequency to s . The solutions to the right of the 
dots do not possess an ergoregion, whereas the solutions to the left of the dots as well as in 
the spirals do possess an ergoregion. 

The onset of ergoregions thus almost always occurs on the main branch of boson star 
solutions, supposed to be classically stable. The presence of ergoregions then tends to 
diminish the physically relevant - since stable - set of solutions. As seen from the figures, 
the diminishment is greater for larger gravitational coupling and larger rotational quantum 
number n, and thus greater angular momentum. However, there always remains a part of 
the main branch of classically stable boson star solutions, not suffering from an ergoregion 
instability. 

Let us now discuss the emergence and structure of ergoregions in more detail for these 
boson stars. For that purpose we exhibit the ergoregions of 1+ boson stars at k = 0.1 and 
2+ boson stars at k = 0.1 and 0.2 in Fig. 12. 
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FIG. 12: Location and shape of the axially symmetric ergoregions of boson stars exhibited via 
cross-sections with the xz-plane: 1 + boson stars at k — 0.1 (upper row) and 2 + boson stars (lower 
row) at k — 0.1 (left) and at k = 0.2 (right). The numbers in the legend indicate the respective 
values of the charge Q. 



Beginning with the 1 + boson stars, we observe, that at the critical frequency (where 
Q s=s 986 for n = 0.1) the condition (46) is satisfied precisely on a circle in the equatorial 
plane. In the (upper left) figure, only the two points of intersection of this circle with the 
zz-plane are seen. As Q is increased (i.e., u s is decreased) the circle turns into a torusdikc 
surface. In the figure, its intersection with the zz-plane is then seen as two small circles 
(Q = 911). From this approximately maximal value of the charge, the set of solutions now 
enters the first part of the spiral, where stability is lost anyway. 

As the solutions evolve into the spiral the charge decreases. But the ergoregion continues 
to grow in size, while it remains (approximately) centered at the location, where it first 
arose. The ergoregion then tends towards its largest size (Q ~ 800), while the center of 
the torus of the ergoregion starts to shift inwards. Interestingly, the ergoregion is not only 
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topologically a torus, but also geometrically almost a torus. 

As seen in the (upper right) figure, with further decreasing charge Q, the ergorcgion then 
starts to decrease in size again, when the solutions evolve deeper into the spiral. At the same 
time, the center of the torus of the ergoregion continues to move further inwards. Thus the 
innermost part of the ergosurface gets increasingly closer to the origin, but never touches it. 

The ergoregions of the set of boson stars at different values of the gravitational coupling 
constant k exhibit an analogous behaviour. Likewise, the ergoregions of the sets of boson 
stars at higher n exhibit an analogous behaviour in the regions studied. As seen in the 
(lower part of the) figure, the shift of the centers of the tori of the ergoregion arises faster 
for n = 2 than for n = 1, yielding almost touching ergosurfaces already in the vicinity of 
the maximal size of the ergosurfaces. 



2. Negative parity boson stars 



For boson stars with negative parity we observe analogous features as for positive parity 
boson stars. This is illustrated in the following. We first consider l~ boson stars for 
gravitational coupling k = 0.1 and 0.2. The dependence of their charge Q and mass M on 
the frequency uo s is exhibited in Fig. 13. For k = 0.1 a large part of the spiral could be 
obtained numerically. 
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FIG. 13: The charge Q (left) and the mass M (right) versus the frequency lo s for 1 boson stars 
for gravitational coupling n = 0.1 and 0.2. The dots indicate the onset of ergoregions. 



The mass M of negative parity boson stars exhibits the same characteristic cusp structure 
as the mass of positive parity boson stars, when considered as a function of the charge Q 
[5]. This is seen in Fig. 14 for 1~ boson stars at gravitational coupling k = 0.1 and 0.2. 




FIG. 14: The mass M versus the charge Q for 1 boson stars for gravitational coupling k = 0.1 
(left) and 0.2 (right), together with the mass for Q free bosons, M = ttibQ- 



In Fig. 15 we compare the 1 boson stars to the 1 + boson stars. In particular, we exhibit 
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the charge Q versus the frequency uj s for gravitational coupling n — 0.1 and 0.2. 
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FIG. 15: The charge Q versus the frequency uj s for 1 boson stars and 1 + boson stars for gravita- 
tional coupling k = 0.1 (left) and 0.2 (right). 

The n-dependence of the negative parity boson stars is demonstrated in Fig. 16. Here 
the charge Q is exhibited versus the frequency u) s for coupling constant k = 0.2. The figure 
also shows the mass M versus the charge Q for 2~ boson stars for k = 0.2, exhibiting part 
of the cusp structure. 




- M.M(Q) 
M = m B Q 




FIG. 16: Left: The charge Q versus the frequency u) s for l - and 2~ boson stars for gravitational 
coupling k = 0.2. The dots indicate the onset of ergoregions. Right: The mass M versus the charge 
Q for 2~ boson stars for gravitational coupling k = 0.2 together with the mass for Q free bosons, 
M = rn B Q- 



Considering the spatial structure of the boson stars with negative parity, the scalar field 
(j) here typically gives rise to a double torusdike energy density T u for the 1~ and 2~ boson 
stars on the main branch, seen already in the negative parity Q-balls. As in flat space, 
the radii of the double tori of the energy density increase when n and thus the angular 
momentum increases. In the spiral the spatial structure of solutions can become more 
complicated, exhibiting for instance a quadruple torusdike structure. 

Addressing the occurrence of ergoregions for the above sets of boson star solutions, we 
observe that rotating boson stars with negative parity also possess ergoregions in a large 
part of their domain of existence. The onset of the ergoregions for the 1~ and 2~ boson 
star solutions is again indicated by dots in the respective figures. Again, if the onset occurs 
on the main branch, then the solutions to the left of the dot and in the spiral do possess an 
ergoregion, while the solutions to the right of the dot do not possess an ergoregion. Likewise, 
if the onset occurs only in the spiral, then the solutions further down the spiral all possess 
ergoregions, while the other solutions do not. We demonstrate the emergence and structure 
of these ergoregions in more detail in Fig. 17, where we exhibit the ergoregions of 1~ boson 
stars at k = 0.1 and 2~ boson stars at n — 0.2. 

Addressing first the 1~ boson stars, we observe, that at a critical value of the frequency 
(where Q ss 540 for k = 0.1) the condition (46) is satisfied precisely on two circles, located 
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FIG. 17: Location and shape of the axially symmetric ergoregions of boson stars exhibited in the 
xz-p\ane: 1" boson stars at k = 0.1 (upper row) and 2~ boson stars (lower row) at k = 0.2. The 
numbers in the legend indicate the respective values of the charge Q. 



symmetrically w.r.t. the equatorial plane. As seen in the (upper left) figure, the circles 
turn into two torus-like surfaces, when Q is decreased in the first part of the spiral. The 
ergoregion thus consists of two disconnected tori, after it first arises. 

Following the set of solutions further into the spiral, the two parts of the ergoregion 
increase in size until they touch at a critical value of the charge (Q w 447 for n = 0.1) 
in the equatorial plane, while their centers shift only little. Beyond this critical value, the 
ergoregion then consists of only one connected piece, a topological torus with an 8-shapcd 
cross-section. 

As seen in the (upper right) figure, the ergoregion then rapidly changes shape and assumes 
a more ellipsoidal cross-section, when the solutions evolve deeper into the spiral. At the 
same time, the center of the torus moves further inwards. Again the innermost part of the 
crgosurface gets increasingly closer to the origin, but never touches it. 

The ergoregions of the set of boson stars at different values of the gravitational coupling 
constant k exhibit an analogous behaviour. Likewise, the ergoregions of boson stars at 
higher n also exhibit an analogous behaviour as seen in the (lower part of the) figure. 



VI. CONCLUSIONS 



We have addressed boson stars and their flat space limit, Q-balls, recalling and extending 
known results for positive parity solutions and presenting new negative parity solutions. Our 
main emphasis has been to study the general pattern displayed by these regular extended 
objects, to determine their domain of existence, and to investigate their physical properties. 

Q-balls and boson stars exist only in a limited frequency range. Whereas both mass 
and charge of Q-balls assume a minimal value at a critical frequency, from where they rise 
monotonically towards both smaller and larger frequencies, boson stars show a different 
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type of behaviour. Their mass and charge tend to zero when the maximal frequency is 
approached, while for smaller values of the frequency the charge and mass of boson stars 
exhibit a spiral-like frequency dependence. The spirals end in limiting solutions with finite 
values of the mass and charge, which depend on the gravitational coupling constant n. 

Q-balls and boson stars possess radial excitations, which so far have only been studied for 
spherically symmetric solutions, and even there a systematic study has not yet been per- 
formed. Preliminary results indicate, that radially excited Q-balls and boson stars possess 
an analogous frequency dependence as the corresponding fundamental solutions. 

Rotating Q-balls and boson stars possess an angular momentum J which is quantized in 
terms of the charge Q, J — nQ, with integer n. The positive parity solutions n + exhibit a 
torus-like energy density, while the negative parity solutions n~ exhibit a double torus-like 
energy density, with the two tori located symmetrically w.r.t. the equatorial plane. Clearly, 
the double torus-like structure of the negative parity solutions is energetically unfavourable 
as compared to the single torus-like structure of the positive parity solutions. 

The classical stability of Q-balls and boson stars can be analyzed according to catastrophe 
theory, implying a change of classical stability at each cusp encountered, when the mass is 
considered as a function of the charge. Both type of solutions then possess classically stable 
branches, representing the physically most relevant sets of solutions, which in the case of 
boson stars have potential applications to astrophysics. 

Rotating objects may possess an ergoregion. But for globally regular objects such as boson 
stars the presence of an ergoregion would imply an instability, associated with superradiant 
scattering [13, 22-24] . Thus rotating boson stars should become unstable, when they develop 
an ergoregion. Therefore the possible presence of ergoregions was put forward by Cardoso 
et al. [13] as an argument to exclude boson stars and various other black hole doubles as 
potential horizonless candidates for compact dark astrophysical objects. 

Our analysis has shown that ergoregions can indeed be present for boson star solutions 
on the classically stable branch. But their presence only diminishes the set of boson star 
solutions with possible physical relevance (where the diminishment is greater for larger grav- 
itational coupling and larger angular momentum), while there always remains a part of the 
branch of classically stable boson star solutions, not suffering from an ergoregion instability. 
Therefore, rotating boson stars cannot yet be excluded as potential candidates for compact 
astrophysical objects. It remains to be seen, whether boson stars with appropriate values 
of the physical parameters to fit observational data will or will not suffer from such an 
ergoregion instability. 
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